Thet a- vacuum and large N limit 
in CP^ -1 a models 



M. Aguado 1 * M. Asorey 2t 
1 Max-Planck- Institut fiir Quantenoptik 
Hans-Kopfermann-Str. 1, D-85748 Garching (Germany) 
2 Departamento de Fi'sica Teorica, Facultad de Ciencias 
Universidad de Zaragoza, E-50009 Zaragoza (Spain) 



Abstract 

The 8 dependence of the vacuum energy density in CP^ -1 models is 
re-analysed in the semiclassical approach, the 1/N expansion and argu- 
ments based on the nodal structure of vacuum wavefunctionals. The 1/N 
expansion is shown not to be in contradiction with instanton physics at 
finite (spacetime) volume V . The interplay of large volume V and large 
N parameter gives rise to two regimes with different 6 dependence, one 
behaving as a dilute instanton gas and the other dominated by the tra- 
ditional large N picture, where instantons reappear as resonances of the 
one-loop effective action, even in the absence of regular instantonic solu- 
tions. The realms of the two regimes are given in terms of the mass gap 
m by m 2 V <C N and m 2 V S> N, respectively. The small volume regime 
m 2 V N is relevant for physical effects associated to the physics of the 
boundary, like the leading role of edge states in the quantum Hall effect, 
which, however, do not play any role in the thermodynamic limit at large 
N. Depending on the order in which the limits N — Y oo and V — > oo are 
taken, two different theories are obtained; this is the hallmark of a phase 
transition at 1/N = 0. 

1 Introduction 

A lack of nonperturbative analytical methods haunts the study of the infrared 
behaviour of confining field theories such as QCD. The main tools used for this 
purpose rely on approximations (e.g., semiclassical, large number of colours), 
and rigorous results are attainable in few corners of parameter space. The 
bordering region between topology and field physics is especially troubling, since 
different methods arise sometimes from apparently incompatible hypotheses and 
physical pictures. 

Two-dimensional CP^" 1 sigma models [2H1I3] are regarded as a convenient 
testing ground to prepare the assault on four-dimensional gauge theories, be- 
cause both kinds of theories share a number of important properties: conformal 
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Figure 1: Structure of the vacuum energy density in the traditional large N pic- 
ture (above) and in the semiclassical picture dominated by instantons (below). 

invariance at the classical level, asymptotic freedom, dynamical mass genera- 
tion, confinement, existence of a topological term 9 and instantons for all values 
of the number of colours N. 

A relevant problem in these theories with topological properties is the 9 
dependence of the vacuum energy density, the quantity that determines the 
phase structure of the theory (for a recent review, see [1]). In particular, the 
fate of the discrete parity symmetry P upon quantisation at the values 9 = 
and 9 = n (the only values for which it is classically conserved) is an issue. 

Perhaps the simplest model in which the subtlety of the ^-dependence of 
the vacuum energy £q is manifest is the quantum rotor [5], i.e., the quantum 
mechanical problem of the dynamics of a charged particle on the circumference 
S 1 enclosing a magnetic flux 9. In the absence of perturbations, the vacuum 
energy is quadratic in 9] periodicity of the physics in 9 — > 9 + 2tt imposes 
that the ground level is twofold degenerate for 9 — tt (i.e., half a flux quantum 
across the region bounded by S 1 ) and there parity is spontaneously broken. 
This, as we will see, mimicks the traditional picture of the large N expansion in 
CP models. However, even slight perturbations compatible with reflection 
symmetry lift the degeneracy of the rotor, by a level repulsion mechanism, 
making the curve £q(9) smooth at 9 = tt. A convenient approximate method 
is that of the dilute instanton gas, where the vacuum is understood in terms 
of tunnelling processes among classical vacua. In the dilute approximation, the 
vacuum energy (a pure nonperturbative effect) is a smooth periodic function of 9 
proportional to (1 — cos 9). This corresponds to the semiclassical approximation 
in CP^ -1 models, where instantons play an all- important role. 

These two regimes have the following paradigmatic expressions for £q(9), 
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illustrated in figure [JJ 

S (9) cx niin{(6» + 27rfc) 2 ; k € Z} (large JV), (1) 

£o(#) oc 1 — cos6> (semiclassical). (2) 

We now consider the situation in CP^ -1 models. Exact solutions are known 
for the quantum CP 1 model (equivalent to the 0(3) model) both at 9 = and 
9 = 7r. In the first case [5], the solution exhibits a mass gap, the spectrum 
consists of an SU(2) triplet, and parity is conserved. This agrees with the 
Haldane map [7], which transforms this model into a chain of integer classical 
spins. Vafa and Witten [S] argued that there is no first order phase transition 
with spontaneous parity breaking at 6 = for QCD, their argument being 
applicable straightforwardly to all CP^ -1 models (see [5] for a proof of the 
Vafa- Witten theorem using the topological charge as an order parameter). 

The exact solution of the quantum CP 1 model at 9 = tt JO] also conserves 
parity but shows no mass gap (this result was anticipated in [TT]). The critical 
behaviour of the model is described by an SU(2) WZNW model at level k = 1. 
This also agrees with the Haldane map, which transforms this model into a 
chain of half-odd spins. By the Lieb-Schulz-Mattis theorem [T2], the absence 
of mass gap implies that P is conserved. On the other hand, the absence of a 
first order phase transition with spontaneous P breakdown at 9 — tt has been 
argued to hold for all CP^ 1 models [T3"| . by analyzing the nodal structure of 
the vacuum in the Hamiltonian formalism [TJ in analogy with QCD [15] . 

For the intermediate region < 9 < tt, analytical techniques are lacking, 
and we must rely on approximations and numerical simulations. We will dis- 
cuss two important approximations, which have been argued to be mutually 
incompatible: the semiclassical method and the 1 /N expansion. 

The semiclassical approach J16 is based, as in the case of the rotor, on the 
picture of the quantum vacuum of Ad gauge theories and 2d CP^" 1 models built 
from tunnelling processes among classical vacua. These nonperturbative pro- 
cesses are dominated by instantons and antiinstantons, (anti)selfdual solutions 
of the classical Euclidean equations of motion. A dilute gas approximation gives 
a 9 dependence of the vacuum energy density of the form 

£q(6) cx m 2 (1- cos 9), (3) 

where m is the mass gap. This dependence cannot be seen in perturbation 
theory due to the nonanalytic dependence of the mass gap on the coupling. 

However, a vacuum based on a dilute gas of instantons and antiinstantons 
is not satisfactory, since the statistical ensemble is dominated by the infrared 
divergent contribution of arbitrarily large instantons, whose density n as a func- 
tion of size p is 

n(p)dp cx (Ap) N ^ (4) 

for the CP^ -1 model, with A a typical scale of the theory. A statistical me- 
chanical treatment of interacting instanton fluids has been developed [T7J [TB] , 
bringing about the instanton liquid picture of the QCD vacuum |19j . This may 
be very relevant for the behaviour of these theories at finite temperature and 
high density. We note in passing that the dilute gas approximation breaks down 
as well in the ultraviolet for the CP 1 model, as pointed out by Liischer [20 build- 
ing on his work with Berg |21) on the geometric definition of a topological charge 
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density on the lattice. Technically, the topological susceptibility in this model 
does not scale according to the perturbative renormalisation group due to small 
distance fluctuations. This is reflected in the singularity of (01 as p — > for 
N = 2; it may also be understood a consequence of the slow vanishing rate of 
the density of Lee- Yang zeros as 8 — > [9]. What is remarkable is that, although 
the semiclassical analysis does not reveal any ultraviolet instanton singularity in 
the case of CP 2 model, numerical simulations suggest a similar pathology [22] . 

The 1/N expansion [23] [23] stands as an alternative to the semiclassical 
method. This technique is based on the simplification of both Ad SU(iV) gauge 
theories and 2c? CP^ 1 models when TV is taken to infinity keeping certain 
parameter combinations constant. 

The 1/N expansion of CP^ 1 models, as developed in [3J and [25], agrees at 
8 = with the known spectrum, given by a massive particle in the adjoint rep- 
resentation of SU(iV). The mass mx is generated dynamically, and the particle 
turns out to be a composite state of two fundamental fields, bound together by 
a Coulomb potential. For 8 ^ 0, this analysis predicts a quadratic 8 dependence 

of the vacuum energy density around 8 = 0. This dependence can be made 
to agree with the fundamental requirement that physics be periodic in 8 with 
period 2tt only if there is a first order cusp at odd multiples of 8 = tt, i.e., a 
first order phase transition accompanied by spontaneous parity breakdown, as 
shown in the upper part of figure [1] 

Instanton effects, being nonpcrturbative, are not visible in the perturbative 
expression ([5]). This led Witten [25] to argue that the 1/N expansion is not 
sensitive to instantons — equivalently, that instantons play no significant role 
in the quantum CP^ -1 models (or in Ad gauge theories) to the extent that the 
1/N expansion is a good approximation thereof. Jevicki [25], however, argued 
that instantons resurface in the 1/N expansion as poles of the integrand of 
the partition function Z (eqn. [5]), and that Z can be computed both by the 
saddle point method and by using a functional Cauchy theorem summing the 
residues of all these poles (representing resonances). Then the large N limit and 
instanton effects would not be a priori incompatible with each other. 

The quadratic dependence §5§ agrees with the holographic picture provided 
by the Maldacena conjecture [27], and moreover with lattice measurements of 
the topological susceptibility of CP N ~ X models (see [4] for a review). The 
Witten- Veneziano formula [28] 129]. derived in this approximation, gives a phe- 
nomenologically correct value of the rf mass in terms of the topological suscep- 
tibility at 8 = 0. However, the appearance of a first order cusp at 8 — tt is in 
contradiction with the results arising from the nodal analysis of the vacuum [13j , 
and with the intuition that level repulsion generically destroys level crossings. 

In this work we show how this discrepancy stems from the fact that the 
large N limit and the thermodynamical limit do not commute. The traditional 
formulation of the 1/N expansion starts directly at infinite spacetime volume 
V = LT = oo. As we shall see, a procedure in which the thermodynamic 
limit is taken after the N — > oo limit provides results compatible both with 
instanton physics and with the rigorous results at 8 — tt, and different from the 
reverse order of limits. A finite volume analysis is in order. This agrees with 
Schwab's [30] |3T] and Minister's [32] [33J approach in the case of the sphere; 
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we moreover outline the application of Jevicki's residue method. We find the 
case of the torus much more tractable and amenable to explicit computation 
after integration over the dual torus parametrising the different holomorphic 
bundle structures within each topological charge sector. In particular, Jevicki's 
approach requires computation of the residues of meromorphic functions instead 
of functionals, and his programme can be carried out in the simplest cases 
exhibiting the reappearance of instantons as resonances (poles) in the one-loop 
effective action. Remarkably, in spite of the absence of regular unit charge 
instantons on the torus |34) . the contribution of this sector is nonzero in the 
resonance approach: this we interpret as the effect of rough configurations near 
the forbidden regular instanton. 

At finite volume there are two regimes: one dominated by instantons for low 
mass theories, m 2 V -C N, and another regime where they are are strongly sup- 
pressed, m 2 V 3> N. The second regime is the relevant one for CP^ -1 theories 
in the thermodynamic limit, but the other regime is relevant for effects where 
the finite volume or space topology play a leading role, like in the appearance 
of edge states in the quantum Hall effect. 

The structure of the article is as follows. In section [21 the traditional large 
N picture of CP^ -1 models is reviewed. The 9 dependence of CP^ -1 models 
formulated on the sphere is considered in section [21 The corresponding analysis 
for the case of the torus is performed in section [4) The consequences of this 
analysis are discussed in section [5j 



2 The traditional picture of the 1/N expansion 

The traditional large N picture of CP^ -1 models was developed in [3] and (25) . 
We shall now give a brief account of it before the analysis in finite volume. 

The large N method is based in a saddle point approximation of the par- 
tition function, defined on the infinite 2d Euclidean plane, after integration of 
the fundamental "J, ^ fields (taking values in C^, i.e., in representatives of 
projective classes in CP^ -1 ). We introduce the dummy U(l) gauge field 

Au = ~ \ (* f M * ~ (^*) f *) . ( 6 ) 

and a scalar field a(x) imposing the constraint = 1 at each point as a 

Lagrange multiplier. Starting from the full partition function at 9 = 0, 

Z = [wV&VAn S[^V - ll exp j- — * [ d 2 a: |D M *| 2 1 

= J V^V^VAf, Va 

xexpj-^/ d 2 x\D^\ 2 -^[ d 2 xaW(ftf-l)}, (7) 
I z 9o Jm 2 z 9o Jm. 2 ) 

we perform the Gaussian integration over "I^ to obtain 

Z= fvA ll Vae- NS ° tt[A '" a] , (8) 



■5 



where the effective action is 

S eS [A^ a] = Tr In (-D* - a(x)) - ^ j_J?xa{x). (!)) 

The saddle point equations 



2^o JR ! 



to(z) -D£ + a 2 ff ^ 

2*^— (*,*) = <), (11) 



a 



can be solved within a renormalisation scheme to yield a saddle configuration 
A fl =0, a = t?4 = M 2 ex P I TTT 1 ' ( 12 ) 



&(/*)/' 

where fi is a mass scale and the corresponding renormalised coupling. 

Perturbation theory around the saddle configuration reveals a dynamical 
system of an iV-plet of charged scalars ^, with a short range interaction due 
to the field a, and electromagnetic interaction due to the field A^. The latter 
develops an effective kinetic term and couples to the scalars with effective electric 
charge e e g = y/Yhvmf^ / N. Thus there is a confining Coulomb interaction 
between scalars, and the spectrum at 6 = consists of (^vj/t) bound states, 
with mass gap 2mx, living in the adjoint representation of SU(iV). 

For 9 ^ 0, the topological term in the action is 

-iOQ = -i— J d 2 xF 01 , (13) 

Q being the magnetic flux associated with the U(l) field A^ and its field strength 
F^, or equivalently, the topological charge of the field ^, which is an integer 
for smooth finite-action configurations. This term plays the role of an external 
electric field in electrodynamics. Therefore, in this picture of the CP^ -1 models, 
its contribution to the vacuum energy density is 

yielding a topological susceptibility 

This quadratic dependence is perturbative, i.e., it can be seen in terms of Feyn- 
man diagrams. Instanton effects, nonperturbative in nature, were argued in |25j 
to be exponentially suppressed in the 1/N expansion, and therefore irrelevant 
for the physics of the CP^ -1 models. However, we have seen that the level 
crossing and first order phase transition at 9 = n implied by (|14[) and the re- 
quirement of 27r-periodicity in 9 are in contradiction with the nodal arguments 
of [13] ■ We will next go over to a compact space with the purpose of showing 
that this incompatibility stems from the infinite volume starting point of the 
traditional 1/N analysis. 
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3 1/iV expansion on S 2 



Clarifying the interplay of N and the volume, and the effects of taking these to 
infinity in different orders, requires the CP^" 1 models to be first formulated in 
a compact (Euclidean) space. 

Schwab (301 ED] and Minister [H|33] studied the 1/N expansion of CP^ 1 
models on S 2 , and observed that the k = 1 contribution to the partition function 
is dominated, for large N, by a saddle point given by a rotationally invariant 
instanton (in the sense that global U(iV) transformations can be compensated 
by 0(3) rotations in Euclidean space). The saddle point equations 

1 ( ] = J_ 
— D^D^ + a 2. 9o 2 ' 

°" ~(x,x) = Q, (16) 



admit for N > |fe|, in a uniform topological charge density background, solutions 
with constant a(x). Indeed, the second equation holds due to parity. The first 
equation states rotation invariance of the propagator G{x, x) of a particle with 
mass \fa. It is easy to show that G(x, y) depends only on the geodesic distance 
between x and y, and therefore the first equation has solutions with constant a. 

3.1 The effective action on the sphere 

Thus, we begin [351 136) with a spherical spacetime of radius R and volume 
V = 4irR 2 , and the action of the CP W_1 model on a background of topological 
charge k, 

S * = ~ T2 I * tA ^ + TT2 / ™ 2 - 1) . ( 17 ) 

^9q J z 9o J 

where integration implies the measure d 2 a; y/g, and is the covariant Laplacian 
in the background chosen. The magnetic flux for the composite U(l) field is 
quantised, 

$ B = 4irR 2 B = 2irk, k e Z. (18) 

We rewrite the constant saddle point value of the a field as to 2 , variable still to 
be integrated upon. 

Integrating out Vt, fy* yields the functional determinant of the operator 
— Afc + to 2 , which is computed in the £ function renormalisation scheme at 
energy scale \x. Discarding unessential factors, 

dm 2 e- NS ^, (19) 

with effective action 

Sf ^ndet^-A^-j-^TO 2 

, , , 2 

= In dct j\. — — — — 2 — to . (20) 
2.9o 
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The eigenvalues of A are 
A n = 



H 2 R 2 

An 



1 n + 



m 2 i? 2 



n = 0, 1, 2 



, _l , . . . , 



(21) 



with degeneracy d n = In + \k\ + 1. 

We use the ^-function definition of the determinant, equivalent to a renor- 
malisation at scale /x: 



In det = rf„ In A„ = rf„ In A„ — I d„ ) h\{fx 2 R 2 ) 
n=o n=0 \n=a / 

^-C~(0)- CiWhVi? 2 ), 



(22) 



being the £ function associated with the operator A = n 2 R 2 A, i.e. the 
analytic continuation of 



n=0 



(23) 



for all complex s ^ 1. 

From the small s expansion of (|23|) . we obtain the effective action for topo- 
logical sector k: 



sf = -CA0)-U0)Mn 2 R 2 ) 



k 2 1 

T + 4 



4vri? 2 
"^1 



m 2 i? 2 + m 7 i? 2 - - ln( M 2 i? 2 ) 



AttR 2 



2\l— + --m 2 R 2 In- 



|fc|+i 

2 



fc 2 | 1 



m 2 R 2 



1*1+1 



A; 2 i 1 



— m 2 R 2 



2Ch -i; 



|fc| + l //c 2 1 



1 m 2 i? 2 

4 4 



(24) 



Here we have used the Hurwitz zeta function £jy(s;u), defined by analytical 
continuation to all s =/= 1 of 



Cff(s;v) = 5Z(n + w) s = 



1 



r( s ) 



oo , s — i v t 



dt 



1 -e- 



-, Res>l, (25) 



and its derivative ^(s;u) with respect to s. Function ([M]) is defined for all 
complex values of m 2 R 2 , bar isolated singularities. 
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Figure 2: Effective action in the k — 1 sector on the sphere, for real m 2 R 2 . 



3.2 Zeros and saddle points of the effective action 

In order to compute Zk, integration over m 2 is still to be performed, through 
imaginary values in order to ensure convergence. 

Let us study the behaviour of the effective action for reai m 2 R 2 (see hgure 
El). To begin with, the integrand of Zk has N(2n + \k\ + l)-fold poles at 

m 2 R 2 = Pn = - (n+^ + lj (n+^j + n = 0, 1, 2, . . . , (26) 

reproducing the eigenvalues and degeneracies of — +m 2 . In this sense, in the 
way predicted by Jevicki, the partition function can be computed by deforming 
the integration curve so as to surround the poles, and summing the residues at 
each of them. However, the problem of computing and summing the residues is 
too difficult to be tackled analytically, although the previous formulae can be 
used in a numerical approach (more progress can be made analytically in the 
case of the torus, as will be seen in next section). 

Alternatively, we can use the saddle point method. The zeros of the deriva- 
tive, 



(27) 

alternate with the poles in the real m 2 R 2 axis, as seen in figure [5] (here i/j(z) = 
T'(z)/T(z) is the digamma function). There is a unique saddle point sq to the 
right of the first pole Po = — ^r, which we assume to be dominant. 

The partition function of sector k in the saddle point approximation is 

k ~ R 2 ]j N\Sf"(s )\' [M) 



dSf 
d{m 2 R 2 ) 



2r>2\ 



Air 



1> 



1*1 + 1 



k 2 + 1 



m 2 R 2 
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up to quadratic order. Explicit results can be obtained for large m 2 R 2 , in which 
region the effective action can be expanded as 



Sf = - ( m 2 R 2 - I) In ™L - ^m 2 + m 2 R 2 
1 3/ a* 2 2^ 

A: 2 1 \ 1 fk 2 4 \ 1 



2A 15jm 2 R 2 \40 315 / m 4 i? 4 
and the saddle point is found to be 

1 fk 2 



O (m- e R- 6 ) , (29) 



m\R 2 = m\R 2 



3 V 24 90 J m\R 2 



where m 2 -* = /i 2 exp {— 2^/^^} is the infinite volume saddle point in (fT2")l . 
The total partition function after summing all topological sectors is 



fee 
fcez 



/2tt m T f 2ttN , , TV 

— — exp <^ — ^- - Nm^R + nn „ n 9 
AT i? I 3fif| 1 gOm 2 ,^ 2 



X 

fee 



2tt m T f 2vrA^ , , TV 

exp <^ — = NmiR 2 H = — - 

TV i? F \ 3g% T OOm 2 ^ 2 



x 3 



2tt 



247rm 



i^)(l + C»K^)), (31) 



where the last equation uses Jacobi's $3 function: 

■& 3 (z\t) = e i,rrn V 27m *. (32) 

Two asymptotic regimes for (f3"Tj) can be analysed. For TV » m^i? 2 , the 
sum therein can be truncated, keeping just the k — —1, 0, 1 sectors. Then the 
vacuum energy density has a typical dilute instanton gas 9 dependence, 

1 z (s «H9) 

ex P I - D a \ i 1 ~ cos )- ( 33 ) 



AirR 2 * \ 24m 2 R 2 



But if m\R 2 3> A^, using the Poisson resummation formula for the 9 function 
in (|31[) and keeping the dominant term in the dual sum, we have 



7 (s o)m „ 4V3^m 2 [2^iV , 2 /V 6m 2 i? 2 ~ 2 | 

Z (0) iV - eXP \^| _JVmTii + 90^p 2 "" ^V~ 6 /' (34) 
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9 being the angle in (— 7r, +tt] differing from 9 by an integer. The corresponding 
vacuum energy density coincides with the traditional large N prediction: 

fW W -fW(0)«M? 2 , (35) 

which is periodic in 9 and undergoes first order phase transitions with level 
crossing at 6 = (21 + 1)tt, £ E Z. 

Before commenting on these two different limiting procedures, let us perform 
the same analysis on the torus [35J [35] ■ 



4 1/N expansion on T 2 

We consider a toric spacetime of linear size L and spacetime volume V = L 2 . 
Functional integration over the fields of the CP^ -1 model on the torus in- 
volves an additional variable, the complex coordinate u <E T 2 in the dual torus 
parametrising the different holomorphic bundle structures associated with the 
complex line bundle E k (T 2 , C) [37]. 

In this case, the effective action Sjf 1 [A^ , a] resulting from integration of the 
ty, ^ fields does not have saddle points. Specifically, the saddle point equations 



(x, x) = 



■2- 



D 2 +« V ' ; 2 5o 2 
D ^ — {x,x) = 0, (36) 



-D2 + o 



do not have solutions with constant a and topological charge density. 



4.1 Quantum saddle points and the effective action on T 2 

The arguments used for the sphere can, nevertheless, be adapted to the torus, 
generalising the saddle point method. By integrating over u (i.e., averaging over 
u), the irregularities of the saddle point configurations are swept off [35l 136] . 
Upon integration over u, a reduced effective action SVed obtains, 

exp{-5 red [^, a, N}}= [ d 2 ue- NS ^^ a \ (37) 

which can be argued to be dominated by constant topological density in the 
large N limit. The generalised saddle point equations 

S d 

SredpV a, N] =0, 



5A^ dN 
5 d 
da dN 



S rcd [A^ a, N] = 0, (38) 



hold in this case because they remain finite as N — > oo. Their solutions we call 
quantum saddle points. 
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To compute the effective action in the sectors of nonzero topological charge, 
the £ function method is used, renormalising at energy scale fj,: 



,,2r2 ™2r2 

--dm-cKfn^.!^.. (39) 

The spectrum of the Laplacian on the torus, in a background with uniformly 
distributed topological charge 2nk ^ 0, is independent of the holonomies u and 
consists of Landau levels — 2oj(n + 1/2), n = 0, 1, . . ., with uj = \B\ = , 
where L is the linear size of the torus. These levels have |fc|-fold degeneracy 
[35]. The zeta function for operator B = (- L 2 A k + m 2 L 2 ) / (47r|fc|) is 



OO 



«')"El*l(" + ^ + 5)"'-l*l<-(^ + i 

yielding the effective action (see [39]) 

ff m 2 L 2 J27T in\k\\ f 1 /m 2 L 2 1 



m 2 L 2 , 47r fc ,,, , f 1 /m 2 L 2 l\l 

= -— i ln ^T¥ - fc !ni^r — — + -H 41 

4tt m 2 ^ 2 \ V 47r l^l 2/ J 

where m 2 is the (constant) saddle point value of the a field, and gn = <7/j(/i) 
is the renormalised coupling at scale /i. In the last equation, = [i 2 exp { — 
27r/gfj.(/i)} stands for the large N dynamically generated mass at infinite vol- 
ume. 

Expression (1411) can be checked to coincide with the dominant term in a 
large volume, constant B expansion of the corresponding effective action (|2"4"]l 
for the sphere: 

S!l k V ^S% k + O(V ) (B= const). (42) 

The contribution of topological sector k to the partition function on the 
torus now depends only on m 2 , all other fields having been integrated out: 

Z k = [ dm 2 e- NS *. (43) 



As in the previous sections, the integration is performed through imaginary 
values of m 2 to guarantee convergence. 

In order to make the functional dependences in some the following expres- 
sions clear, it is useful to define dimensionless variables 

m 2 L 2 m^L 2 

v= Mky y ° =4^r (44) 
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S(y) 




y 



Figure 3: Effective action S(y) for the torus, with y 6 K. 



Then the fc-sector partition function is 
/ i i\\ JV l fe l 

Z k = ^fdy e--l^-- = ^ / dye-WM, (45) 

where the function 

%) = ^=^ -ln r( ^ /2) (46) 

|«| V27T 

is defined for all complex values of y, except for a series of poles of the integrand 
of Zk , as can be seen in figure [31 

For small values of y, that is, when \k\ <C m 2 L 2 , the exponent simplifies: 

§ ^ = -y^i + y + iy-2sh +0 ^ 4 ^ ^ 

meaning that the effective action has an expansion in powers of the topological 
number k where the first nontrivial term is quadratic: 

ff _ m 2 L 2 m 2 m 2 L 2 irk 2 f fc 4 \ 

In the opposite limit, when |fc| 3> m 2 L 2 , that is, for large y, 



In 2 , , , ,.. 7T 2 „ 1 



- ^ + {lnj/ -V(l/2)}t/- ^y 2 - ^"(l/2)j, 3 - ^2/ 4 + 0(2/ 5 )- (49) 

Written in terms of the effective action, we see that the leading term is linear 
in the absolute value of the topological number: 

off In 2, , m 2 L 2 f m 2 L 2 , , 1 fm 4 L 4 \ , s 

Notice the change of asymptotic behaviour of the effective action S% s in the 
different topological sectors. According to equations (|4*5| and (|50)) . for small 
values of the topological charge, < m 2 V, the effective action is quadratic in 
k, whereas for large topological charges its leading term is linear in \k\ 9- This 
change of asymptotic behaviour has important physical consequences. 
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■S3 



S2 



Si 



■so 



P3 = I{Uf 
P2 = I{II) 2 
Pi = I(II) 

Po = I 



Figure 4: Structure of poles p n of the effective action (squares) and saddle 
points s n of the integrand of the partition function (circles) on the torus, in the 
complex plane of the variable y — ni^L 2 / (4ir\k\). Note that zeros and poles 
tend to coalesce as n grows. 



4.2 Zeros and saddle points on the torus 

Extrema of S(y) allow us to perform a saddle point approximation by deforming 
the integration contour so that it passes through the dominant extremum. The 
exponent in the Zk integral has an overall N \k\ factor, therefore the saddle point 
approximation is a large iV|fc| expansion, and results in terms of y are general 
for all k 7^ sectors. 

Poles of the integrand of Zk give us a chance of testing Jevicki's proposal, 
since functional integration has been reduced to integration along a path in the 
complex plane. The integration contour must be deformed so that it surrounds 
each pole. In spite of the fact that various fields have been integrated out in 
the effective action we are working with, we shall see that instantons reappear 
in these poles. 

The structure of saddle points of Sk and poles of the integrand of Z^ is 
represented in figure [4] 



4.2.1 Poles and Jevicki's approach 

Let us consider the poles. There is an infinite series of iV|fc|-fold poles at values 
m n of m such that 

V=^=** = -(»+\)> " = 0,1,2,.... (51) 

Their multiplicity is equal to the complex dimension of the moduli space of 
charge fc instantons in the CP^ -1 model (see |40) ) . 
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The partition function can be written as a sum over residues a la Jevicki, 

JV|fc| 



L 2 ^v^p n \ J2n 



e -N\k\ylny 



n=0 

JVIfel oo 



^ (£) 3 E^ |fe|1 - ^ {^r(-» + e )}"'\ (52) 



n=0 



There is no difficulty in computing and summing the residues for the first cases, 
iV|fej = 1 and iVjfcj = 2: 

/ — 27TTOT f mLh 2 1 , , 

^1*1=1 = V2-^exp --J r (53) 



and 



2miL 2 



^| fc |=2 = m|K (^j_j (54) 

(where Ko is a modified Bcsscl function), but the partition function for higher 
values of iV|fc| turns out to be more difficult to compute. From the expansion 
of (l52"|) . we can write it as 

Zk = —{v^ M~^J T m-iM- (55) 

The function 7jj jTI (j/o) is given by 

OO 

7fl,„(t/o) = <Ws+t,-R a r fo s ct, n , (56) 

r,s,t=0 

with coefficients defined by the expansions 

oo 

Vq 6 = ^2 a r£ r , 



£ 6 - 



' X use- j 



Sin 7T£ 

s=0 



r(» + i-e) ,_„ 

Expressions for a r , 6 S are readily found, 

(-my r 



2(2 S -! - lW'lBj 
: , s even, 



0, s odd, 



(58) 



where B s are Bernoulli numbers. As for C{, n , it can be written as a sum over 
Young tableaux of order t, 

£( - 1 ).-E;.,,nU[f-'(n + D]" 
Y.T.(t) 11^=1 *- ^- 
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where i>j, j — 1, . . . , t, are the numbers of rows with j elements, such that 
y~lj—i jfj — t, and ip is the digamma function. 

As an example, for the k = 1 sector in the CP 1 model we need 

T hn = V>(n + 1) - lny 0j (60) 

from which equation ([54"| obtains. For the k — 2 sector in the same model, 
expression (|59|) is already too cumbersome to compute (|56l) explicitly: 



Ta, n = - ^ (Inyo) 3 - y hryo 



^(n + 1) 



(Inyo) Qv(^ + l) 2 -^'(« + l)) 



+ i V(n + I) 3 - o + 1)VA" + 1) + 7 ^"(n + 1). (61) 
6 2 6 

However, it is not necessary to perform the summation in (|55p to realise that 
the pole structure has a natural interpretation in terms of instantons. From the 
original classical action 

S[tt, tft, ^ a] = ^ / d^lD^p + Z. /" d 2 xa(x)(*tvl/_i), ( 6 2) 
^ffo Jt 2 z 9o Jt 2 

the classical equations of motion 

(-D 2 + a) * = 0, a = I^D 2 * (63) 

ensure that, for classical solutions, the value of the action is given by the integral 
of —a: 

N f N f 

S <* = oH dx \ D ^\=^ d 2 x(-a). (64) 
z 9o Jt 2 z 9o -It 2 

The nth pole of the integrand of Zk corresponds to a value of a — m? such that 

m 2 T 2 f 1 \ N N 

" " =-[n+-) =>■ 5 c i = -^m 2 L 2 = (l + 2n)-^2^|fc|, (65) 



47r|fc| V 2; 2y 2 v '2y 2 

i.e., it exactly matches the classical action of a multiinstanton configuration 
composed of an instanton and n instanton-antiinstanton pairs (figure U) . This 
is in contrast with the case of the sphere, where the structure of poles of the 
integrand of Zk does not correspond to charge k multiinstanton configurations. 

Notice that, although there are no unit charge instantons on the torus [34] . 
there is a nontrivial contribution of the k = ±1 sectors to the total partition 
function. This is so because the dominant configurations in the partition func- 
tion are not the smooth classical solutions, but rather the rough configurations 
in the neighbourhood of these. The classical solution fails to exist for unit 
charge, but this is a classical accident, which does not affect the quantum dy- 
namics of the system because of the nonzero contribution of the neighbouring 
configurations. 
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4.2.2 The saddle point method 

Now let us consider the extrema of S(y). These are the zeros of its derivative, 

dS(y) 



= Inyo ~i>{y + 1/2), (66) 



and constitute a sequence y = s n , n = 0, 1, . . . Saddle points and poles alter- 
nate, so > po > si > pi > ■ ■ ■ , as seen in figures [3] and S) The saddle points 
approach the poles for large n, lim„_ i . 00 s n /p n = 1- 

If the dominant saddle point so lies in the region !/ > 1, we can find its 
location as an expansion in powers of y^ 1 starting from (|47[) : 

1 1 9f)Q7 

*o = l/0 + ^-- i^P- 3 +0(y^). (67) 
24 yo 576 

Equivalently, the infinite volume value m 2 -. of the saddle point receives finite 
volume corrections, 

2 7 r 2 |fc 

We evaluate the partition function in sector fc, up to quadratic order, 

{so) 4^ r*r NW s (so) ( f§\ 1/2 m 

k L 2 y iV|fc| ^dy 2 ^ ' [ ' 

from which 



L 2 y JV|fc| ""^ 1. i ' l '" lau 24?/o ' 16z/o ' 5760yg 

4nm T [ mlL 2 tt Nk 2 ( n\k\ V , , N 1 , . 
"P^-^-fe— fi^+UP 5 ) ) ' (70) 



The partition function in the trivial topological sector can be computed in 
the large m 2 L 2 limit by substituting an integral for the sum in the definition of 
the £ function, 

W^M"^ 1 "^- 1 )}' (71) 

and coincides with the B — > limit of the result for Z^. There are no poles in this 
sector, which is compatible with the absence of instantons of zero topological 
charge, and its only saddle point yields an approximation 

Z^v^exJ-N^), (72) 
V2NL I 4tt J ' 

agreeing with the result ([70|) for the nontrivial sectors. 
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The full partition function in terms of the vacuum angle 9, in the saddle 
point approximation, is the sum 



Z (so) {6) =J2 z i e 
AnrriT 



(s ) -ike 



V2NL 
V2NL 



exp 



exp 



An 

An 



7T Nk 2 

Z_/ ex P 1 g rri^L 2 



E 



iN 



(73) 



Still in the region y ^> 1, i.e. m 2 L 2 3> 47r|/c|, let us analyse this partition 
function in two different regimes, depending on the relation between the number 
N and m\L 2 . 

For TV 3> m^-L 2 , the contribution of high topological sectors (large \k\) to 
(173")) can be neglected. Keeping just sectors = 0, ±1, we obtain 



Attuit 
V2NL 



exp <^ -N- 



r m%L 2 



An 



1 + 2 exp 



7T TV 



6 m^L 2 



cos 6* 



(74) 



giving rise to a vacuum energy density 

£ { Q So) {9)-£ { s °\o) = - 
A 



_L_ Z^(0) 
L 2 n Z( s o)(0) 



L 2 



exp 



7T jV 



6 m^L 2 



(1 - cos 9), 



(75) 



i.e., the typical £q(9) dependence of a dilute instanton gas. This is a 2-7r-periodic 
function, smooth for all values of 9 including 6 = ±tt. Hence, there is no first 
order phase transition. 

This regime would be compatible with a definition of the 1/N expansion 
in which both the large TV limit and the 1/iV corrections are studied in finite 
volume. 

However, if we move to the region m^L 2 ^> N, it proves convenient to use 
the Poisson formula in (1731 to get 



Z {so) {9) = 



Ay/3 



nmr. 



N 



exp <^ —N 



m^L 2 
An 

mLL 2 



i) 



2tt 





N 



f A m^L 2 l x - 
^ ' <?ez 



idtm^L 2 

~~ /v 

Qnm^L 2 
N 



'In 



where the •& function with characteristics is given by 







(z\t) 



exp {inT(n + a) 2 + %2n(n + a)(z + &)} . 



(76) 



(77) 



Now, it is the dual sum in q that is dominated by the low |g| terms. Defining 
again 9 as the angle in [—w, +n] differing from 9 by an integer, 



'( s o) 



(0) 



4%/3- 



N 



exp 



-N- T ! 



An 



2nN 



(78) 
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This reproduces the traditional large N picture, where instanton effects are 
suppressed, and the vacuum energy density depends quadratically on 9 within 
the interval [— w, +%]: 

^) W -fW(0)«M? 2 . (79) 

Periodicity in 9 is guaranteed because £q(9) depends on the periodic variable 
9, but this function is not smooth at odd multiples of n, where 9 is doubly 
defined, levels cross and a first order phase transition occurs. 

This regime is compatible with a definition of the 1 /N expansion in which the 
thermodynamic limit is performed first, and N is taken to infinity afterwards. 

These results agree with the analysis of CP^ -1 models on the sphere. Since 
the physical pictures pertaining to the regimes N 3> m\L 2 and m\L 2 ^> N are 
different, the limits N — > oo and V — > oo do not commute, and the orders in 
which these limits are taken determine different theories. This behaviour points 
towards the existence of a phase transition in the N — > oo theory. 



5 Discussion 



After the careful analysis of the large N method on the sphere and the torus, we 
conclude that the apparent incompatibility between instanton physics and the 
1/N expansion has its cause in the formulation of the latter in infinite volume 
and is a subtle effect of the noncommutativity of the large N and thermodynamic 
limits. 

To clarify this, consider the essential dependence of the vacuum energy den- 
sity on the angle 9, the volume V, and the number of colours N, at fixed saddle 
point mass mx 



£o(0) 



1 

V 



111 



03 ( 


' 6 


iN 




6m^V 


^3 


(° 


iN 


6rr4V ; 



(80) 



which is valid in the cases of the sphere and the torus provided we define the 
free energy as a function of 9 by subtraction of the contribution at 9 = for 
each N. 

The various limits of (180[) are best discussed in terms of a dimensionless 
variable x = N / (Qm\V) and the function 



N 
6m^ 



£ (6) = f(0, x) 



x In 



03 (0 |*as) 



■d 



■ x In • 



0/(2*) 




(oii) 



03 i 



(81) 



The last equation is obtained by applying the Poisson resummation formula, 
that is, the modular transformation of the theta functions. For x strictly posi- 
tive, / is a well defined real analytic function of 9 £ R. The nonanalyticities for 
complex 9 are branch cuts located at the zeros of the Jacobi $3 function, that 
is, for 9 = (1 + ix)n and its translations by integer multiples of 2ir and of i2itx. 
These zeros never occur for real values of 9. 

The two limiting regimes we have been discussing are given in terms of the 
dimensionless quantities as: 
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• Semiclassical: send x to infinity (meaning TV ^> m\V). In this case 

ixj = 1 + 2 e" ra cos 9 + 0(e~ 27TX ), (82) 

and we recover the nonperturbative result 

f(6, x) = xe- nx {l - coa6) + 0(e- 27rx ), (83) 

equivalent to 

1 f TT^V I , ^ /m£ f ttTV 1\ 

=7 exp r 6^7; (1 " cos0) + ° l^r exp r ^)) ' 

(84) 

equivalent to ([5]) (and vanishing as x — > oo together with all of its deriva- 
tives) . 

• Traditional large TV: send x to zero (meaning TV <C m\V). This is a 
problematic limit for the modular parameter r — ix, which leaves the 
upper half plane. We consider separately the regions 9 £ (— tt, tt) and 
9 = tt (the rest of the function is obtained by periodicity): 

/( , 6( ^,..)-£ +0 (.*{-i<^}), 

/(0 = tt, i) = r- -xln2 + 0(e- 7r / x ), (85) 

47T 

reproducing ([1}, including the first order cusp at 9 = tt: 
£ o (0e(-7T, TT)) 
£o(0 = tt) 




^ TO T ^2 

2ttTV 



(9 
In 2 



TV 



exp 



67rm^(l - 9/tt) 2 
TV 



2vrTV n ~ ~V 



O 



TV 



exp 



TV 



(86) 



The role of modular parameter suggests an analogy with finite tem- 

perature models, where the number of colours corresponds to the inverse tem- 
perature p. In this sense, we expect that the physics at 1/TV = corresponds 
to zero temperature phenomena. The fact that the thermodynamic and large 
TV limits do not commute (reflected in the behaviour of x) would suggest the 
presence of a phase transition exactly at zero temperature, i.e., 1/TV = 0. 

Returning to physical quantities, we have shown that finite volume effects 
in the 9 dependence of the vacuum energy density for CP^ -1 models on S 2 
and T 2 , when all topological sectors are taken into account, give rise to two 
asymptotic regimes, one dominated by instanton effects (when TV 3> m^V) and 
the other by the conventional large TV picture (when TV <C m\V). These are 
smoothly connected by an interpolating region. 

It should be realised that the basic hypotheses of the method of large TV do 
not hold when TV <C m^V, for which precisely the traditional large TV results 
obtain. The saddle point technique needs TV to be the largest dimensionless 
parameter of the theory, in particular larger than rrS^V . Two very different 
theories are defined by interchanging the noncommuting limits TV — > oo and V —¥ 
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oo. In principle, the only procedure consistent with the saddle point method is 
taking the large N limit first, and then going over to the thermodynamic limit. 

However, for small values of 9, where the large N approximation is expected 
to hold [H] , both procedures seem to make sense. Lattice measurements of the 
9 = topological susceptibility agree with the traditional large N picture, cor- 
responding to performing first the V — > oo limit, and then taking N to infinity. 
However, these simulations have been only carried out for values of the parame- 
ters such that m\V > N (see Tableland references [H>] [50] [55 ) , agreeing with 
our analysis in the region not validated by the saddle point method. It would 
be interesting to rerun these simulations for smaller volumes, still close to the 
thermodynamic limit with a stable mass gap, but where an instanton-dominated 
9 dependence of the vacuum energy density could emerge and eventually take 
over. Notice that the singularity of the topological susceptibility pointed out by 
lattice simulations in the CP 1 model is also found in the semiclassical scenario. 
The existence of such a singularity can also be understood by the presence of a 
family of Lee- Yang zeros of the analytic continuation of the partition function 
in the complex 0-plane, converging to 9 = in the thermodynamic limit [5]. 



Reference 


771 


L 


m 2 L 2 


N 


m 


Blatter et al. 


0.165 


40 


43.6 


2 


m 


Ahmad et al. 


0.179 


50 


80.1 


2 


E2 


Lian-Thacker 


0.111 


100 


123.2 


2 


m 


Keith- Hynes-Thacker 


0.179 


50 


80.1 


2 


m\ 


Lian-Thacker 


0.084 


100 


70.6 


3 


m 


Campostrini et al. 


0.066 


120 


62.7 


4 


m 


Burkhalter et al. 


0.131 


32 


17.6 


4 


E2 


Lian-Thacker 


0.088 


100 


77.4 


4 


M 


Keith- Hynes-Thacker 


0.180 


50 


81.0 


4 


m 


Ahmad et al. 


0.186 


50 


86.5 


6 


m 


Lian-Thacker 


0.085 


100 


72.2 


6 




Campostrini et al. 


0.196 


72 


199.1 


10 


M 


Del Debbio et al. 


0.192 


60 


132.7 


10 


[43] 


Ahmad et al. 


0.212 


50 


112.4 


10 


m\ 


Lian-Thacker 


0.058 


100 


33.6 


10 


m 


Keith- Hynes-Thacker 


0.212 


50 


112.4 


10 


m 


Del Debbio et al. 


0.397 


42 


278.0 


15 


\M 


Del Debbio et al. 


0.418 


30 


157.3 


21 


M 


Vicari 


0.287 


60 


296.5 


21 


m 


Del Debbio et al. 


0.305 


56 


291,8 


21 


m 


Vicari 


0.411 


42 


298.0 


41 



Table 1: Numerical values of the mass gap m for different CP W-1 models at 
different volumes m 2 L 2 . The models are in the thermodynamic regime m 2 L 2 > 
N in all cases. The role of instantons is only manifest in the cases CP 1 and CP 2 

As regards the neighbourhood of 9 = n, the nodal analysis of [13] appears 
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incompatible with a first order phase transition at 6 — n. This is the behaviour 
of the system for lower values of N for any volume, i.e. CP 1 and CP 2 models. 
For larger values of N, for instance N > 4, this behaviour is only observed for 
small volumes, i.e. volumes which verify m 2 V <C N. For larger volumes the 
effect is swept off by the infrared fluctuations and the system undergoes a phase 
transition at 9 = 7r with spontaneous CP symmetry breaking. 

The behaviour for intermediate values of 9 has so far proved elusive to nu- 
merical techniques, due to the inaccuracies inherent to lattice simulations in 
this region. Let us however remark that a novel technique |52| based on analyt- 
ically continuing the 9 dependence to imaginary values has been introduced to 
overcome this problem. This technique has been applied to the CP 9 model [53] , 
with conclusions agreeing with the usual large N expansion. The consistency of 
this technique for any value of N is supported by the absence of singularities in 
the analytic extension of the partition function to the whole 0-plane [5]. 

Summing up, the large N method is compatible with instanton effects. Be- 
sides, the results on 9 dependence obtained with this tool agree with the known 
behaviour at 9 = 0, i.e., the vacuum energy density is differentiable there and 
the Vafa-Witten theorem holds. It is also compatible with the numerical deter- 
mination of the topological susceptibility at 9 = 0. The analysis of the poles 
of the partition function on the torus supports the method of Jevicki, whereby 
instantonic effects appear in the large N limit in the form of resonances. First 
order phase transitions with spontaneous parity breaking at 9 = ir appear in 
the formulation of the models directly at infinite volume, and we have exposed 
the analytic roots of this fact. The large N method is thus in agreement with 
all exact results on 9 dependence, and provides a valuable bridge between the 
angles 9 — and 9 = tt. 

Let us remark that the behaviour of the theory in finite volume plays a 
fundamental role in condensed matter settings, where sigma models can be 
used as effective theories for the quantum Hall effect [5H [551 [SHI [57] . In this 
context, the Hall conductivity is identified with the coupling of the topological 
term, and the stability of Hall plateaux is linked to the renormalisation group 
running of the couplings (including 9). The large N limit of CP^ -1 models 
was studied in connection with this phenomenon in a series of papers (see, e.g., 
[551 15*51 1501 15T1 IE2]), in which the different regimes we have discussed were also 
identified; in this case the traditional large N limit at infinite volume is blind 
to edge effects, which of course are crucial for the physics of the Hall effect. 
In particular, edge currents are a finite size effect and this suggests that the 
m 2 V <C N regime of CP^" 1 sigma models is the relevant regime for their 
description. 

Finally, we remark that the difference between the two regimes is due to the 
asymptotic behaviour of the effective action in the different topological sectors. 
According to equations (|4"5[) and ([50]) . for small values of the topological charge, 
\k\ < m 2 V, the effective action is quadratic in k, whereas for large topological 
charges its leading term is linear in |fc| [9 . The two regimes also differ at 
finite temperature. Since the spacetime volume is V — LT, the change of 
asymptotic dependence of the effective action on the topological charge can be 
associated with a finite temperature crossover from the low temperature regime 
f3 = 1/T > m 2 L/\q\ to the high temperature regime (3 < m 2 L/\q\ and cannot be 
related to any phase transition [63] . One might expect a similar phenomenon in 
QCD, although in that case there is a finite temperature phase transition [64] . 
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